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Photosynthetic antenna complexes are responsible for absorbing energy from sunlight and trans- 
mitting it to remote locations where it can be stored. Recent experiments have found that this 
process involves long-lived quantum coherence between pigment molecules, called chromophores, 
which make up these complexes Expected to decay within 100 fs at room temperature, these coher- 
ences were instead found to persist for picosecond time scales, despite having no apparent isolation 
from the thermal environment of the cell. This paper derives a quantum master equation which 
describes the coherent evolution of a system in strong contact with a thermal environment. Con- 
ditions necessary for long coherence lifetimes are identified, and the role of coherence in efficient 
energy transport is illuminated. Static spectra and exciton transfer rates for the PE545 complex of 
the cryptophyte algae CS24 are calculated and shown to have good agreement with experiment. 
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When a solar photon is captured by a photosynthetic 
organism, it creates a localized excited state, or exciton, 
in one of several optically active pigment molecules, or 
chromophores, in a specialized antenna complex. The 
exciton is then conveyed via electrostatic coupling to a 
remote reaction center where its energy can be stored 
and ultimately used to do work. Because the exciton 
has a natural lifetime on the order of nanoseconds, it is 
necessary that the transfer process be very rapid, so that 
the exciton's energy can be stored before it decays. 

In order to understand the dynamics of photosynthesis, 
it is necessary to understand the dynamics of the exci- 
ton transfer process. However, such an understanding is 
complicated by the many internal degrees of freedom in 
the chromophores, and by interactions between the chro- 
mophore and the cellular environment. In the limit of 
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FIG. 1: Spatial arrangement of chromophores and protein 
backbone in PE545. 
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weak coupling between chromophores, Forster theory [6| 
describes the transfer of excitons as a fully incoherent 
process. In the opposite limit, of weak coupling between 
chromophores and the surrounding environment, Red- 
field theory [7[ describes the coherent transport of exci- 
tons using a perturbative expansion over a weak system- 
reservoir coupling, but may give unphysical negative or 
diverging exciton populations. Such unphysical behavior 
may be avoided through use of a secular approximation 
[§[ , at the cost of decoupling the evolution of populations 
and coherences. Neither the Forster nor the Redfield 
limits is applicable to the limit in which the coupling be- 
tween chromophores and between chromophores and the 
environment are both strong. 

Nonperturbative approaches to the problem of exciton 
dynamics are less easily classified. The Haken, Reineker, 
Strobl model H 03 calculates the evolution of an elec- 
tronic reduced density matrix due to electrostatic cou- 
plings between chromophores, with dephasing due to in- 
teraction with the environment. This model predicts long 
lived coherences between chromophores, but yields an 
equilibrium state with equal occupation probabilities for 
all chromophores, which is true only in the high tem- 
perature limit. Other approaches include the stochas- 
tic Schrodinger equation [111 [l2[ and methods based on 
quantum walks [13]. Several studies have found that in- 
terplay between the system and the environment, includ- 
ing nonmarkovian effects and details of the system's vi- 
brational structure, can affect coherence lifetimes and the 
rates of exciton transfer [3, [H| ■ 

The need for an improved treatment of exciton dy- 
namics has been highlighted by the recent observation 
of long lived coherences between chromophores in pho- 
tosynthetic antenna complexes [1-4]. Whereas a sim- 
ple energy-time uncertainty approach would suggest that 
such coherences must die in less than 100 fs @, they were 
instead found to persist for picosecond timescales - on the 
order of the time needed for excitons to travel between 
chromophores. Rather than being confined to the low 
temperature limit, such lifetimes were observed even at 
room temperature, shedding doubt on the common as- 
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sumption that biological systems are too "hot and wet" 
to display quantum coherent behavior (l6| . 

The observation of long lived coherences in biological 
systems at room temperature raises several questions. 
First, how can these coherences persist in the thermal 
environment of the cell? Are the coherences somehow 
protected from interaction with the outside world, or 
are they preserved through another mechanism? Second, 
can the organism exploit coherence to improve the effi- 
ciency of photosynthesis? Does coherent transport alter 
the rate at which excitons travel between chromophores? 
Does the quantum information contained in the coher- 
ence terms affect the flow of population through the com- 
plex? Some have speculated that the complex may per- 
form quantum information processing to optimize the 
transfer path although others [ItJ have found 
that such speedup may be short lived. In addition to 
these questions, the underlying problem of a few state 
system with internal degrees of freedom interacting with 
a thermal reservoir may have applications in fields such 
as quantum computing, where loss of coherence due to 
interactions with the surrounding environment is a major 
limiting factor |I8j . 

This paper presents a new theory of coherent evolu- 
tion in a multichromophore system which does not rely 
on a perturbative parameter. Instead, an assumption of 
rapid thcrmalization of vibrational degrees of freedom is 
used to derive equations of motion for a slowly evolving 
electronic reduced density matrix. In this way, both the 
transfer of exciton population and the evolution of co- 
herences are treated on an equal footing. The resulting 
theory gives a simple explanation for long coherence life- 
times in the high temperature limit for a system which 
interacts strongly with its environment. 

Section Q] derives equations of motion for an elec- 
tronic reduced density matrix. Part II CI reduces the full 
electronic + vibrational+environment density matrix to 
yield equations of motion for a reduced electronic density 
matrix, while part II Dl reduces these equations of motion 
to the form of the Haken Reineker Strobl model [H, in 
which both the density matrix and the effective coupling 
incorporate thermodynamic quantities. Section |TT] ap- 
plies the resulting theory to a two chromophore system, 
to find lifetimes for the decay of coherence or the trans- 
fer of population between two chromophores. Section 
IIIII applies this theory to the PE545 antenna complex of 
cryptophyte algae Rhodomonas CS24, shown in Figure 
[TJ comparing to experimentally measured transfer times, 
as well as absorption, circular dichroism and fluorescence 
spectra. 



I. EQUATIONS OF MOTION 

Although the primary function of an antenna complex 
is to transfer electronic excitation between chromophores, 
it is not a purely electronic system in the sense of hav- 
ing electronic degrees of freedom which can be separated 



from the other degrees of freedom in the system. Rather, 
both the ground and excited " states" of a particular chro- 
mophore correspond to an infinite number of vibronic 
states. Each of these vibronic states may interact with 
the local environment of the cell, a system which is com- 
plicated and difficult to characterize. In order to calcu- 
late equations of motion for the purely electronic quan- 
tity of interest, it is necessary to remove the vibrational 
and environmental information from the system - ie, to 
reduce the full electronic + vibrational + environmen- 
tal density matrix to a purely electronic reduced density 
matrix. 



A. The Vibrational/Electronic System 

The full electronic + vibrational + environmental 
Hamiltonian may be partitioned into an electronic + 
vibrational system interacting with a reservoir, H = 
Hs + Hr + HgR. Here the part of the total Hamilto- 
nian dealing with the system is relatively well charac- 
terized, while that dealing with the reservoir is less so. 
To aid in the reduction to a purely electronic density 
matrix, it is useful to write vibronic states as the prod- 
uct of vibrational and electronic eigenstates. Writing the 
system Hamiltonian Hs = Hq 1 + Hy lh + V, system eigen- 
states in the single exciton manifold can be written \i, n), 
where i is the excited chromophore and n the vector of 

excited vibrational states, ft = m, ri2, TI3 In this 

basis, HQ l \i,ri) = Ei\i,n), HQ lh \i,rii) = e,^ \i,rii) and 
e a = E j e «] ■ Vibrational states of the excited chro- 
mophore \i,ni) may differ from those of the unexcited 
chromophore |0, rii) due to differing potential energy sur- 
faces. Here it is convenient to introduce electronic = 
\i) (0| and vibrational 1/+ = \i 7 n>i) (i, 0j| + |0, rii) (0,0i| ex- 
citation operators, so that the state of the complex can 
be written as the product of excitation operators acting 
on an initial ground state 

i,n j 

Deexcitation operators fii and v ni are given by the Hermi- 
tian conjugates of the excitation operators. In the limit 
that electronically exciting one chromophore does not af- 
fect the potential energy surfaces on other chromophores, 
the electronic excitation/deexcitation operators com- 
mute with the vibrational excitation/deexcitation oper- 
ators for a different chromophore. 

B. The reservoir 

In contrast to the system being studied, whose Hamil- 
tonian is relatively well known and studied, the reser- 
voir may be complicated and difficult to characterize. In 
Forster theory, rapid dephasing of electronic and vibra- 
tional degrees of freedom allow explicit consideration of 
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the reservoir to be avoided, as dynamics can be calculated 
from the overlap of the donor emission and the acceptor 
absorption spectra [l9| . In Redfield theory, the reservoir 
is often modeled as an infinite set of harmonic oscillators 




In the current work, each chromophore is assumed to 
interact with a separate reservoir, so that H$r can be 

written as Hsr = Vcr \i) and the system-reservoir 
interaction is incapable of transferring excitation between 
chromophores. As in Forster theory, an assumption of 
rapid thermalization allows electronic equations of mo- 
tion to be derived without explicit reference to the form 
of 

If eigenstates of the reservoir Hamiltonian corre- 
sponding to chromophore i are given by \Q) = 
then the mixing of eigenstates \i,n, Q) is gov- 
erned by V and Hsr. The assumption of rapid ther- 
malization is that the short time evolution of the 
system is given by U(t,t ) w Yl i U io)> where 

U®(t,to) = exp[-i^V«(r)dr]|t){<|. If = 

Si,a,Ci a <(*) 6 ff (*) c Ci(*) N>«)Ci)) the electronically diago- 
nal character of U means that a,(£) does not evolve, while 
the vibrational and environmental degrees of freedom are 
driven to thermal equilibrium 

Tr c<=Ct; 6WC« C «c*P=5 n , n ,pW. (2) 



Here Pk is the probability of vibrational state n on chro- 
mophore i being excited, while — fjj PnJ is the 
probability of a particular vector of vibrational states 
being excited. 

Eq. [2] can now be used to reduce equations of motion 
for the full electronic + vibrational+environment system 
to a purely electronic subsystem. 



C. Density Matrix Reduction 

Electronic equations of motion can be found by tracing 
over environmental and vibrational degrees of freedom. 
To this end, let 

PmCjmr, = «» 6 « C Ci ^> Ci) 0') T]j \ aft^C*. X 

e -i(E t +4+E u -E J ^e ( i i ) -E v .)t ^ 

be the full density matrix written with factorized co- 
effients and 



'WiSA* = IC1) (CI (II b n' 3 K> K \)Pmarnv(R & S K> KIK{ Wf) I ( 4 ) 
n'. m'. 



be an auxiliary matrix defined in terms of the squares 
of vibrational and environmental coefficients, so that 
a i a ] = J2n,mXi,v 3 a inCd"ivy Equations of motion for the 
electronic coefficients a^a* can now be found by substi- 
tuting Eq. [2] into equations of motion for a. 

The rate of dephasing between different chromophores 
is found by substituting Eq. [5] into Eq. |H yielding 

ai a*(t) = a,a*(0)gP| i) P«e- i ^+4 (5 -^-^)*, (5) 

Similarly, equations of motion for coherent evolution 
are found by substituting Eq. [2] into the operator equa- 
tion of motion for a due to Hg + Hr, J^er = —i[Hs + 
Hr,o], yielding 

— aiCij = -i[VifikpPf; 'P^'a k aje y lT « * 1 - 

a^P^^e'^r-^-^ 

(6) 

where Vkpjm = (0| [liV^Vvtp^ |0). 



The complicated forms of Eqs. [5] and [6] which result 
in part from their generality, may be simplified by mak- 
ing assumptions about the vibrational spectrum and the 
form of the vibrational/electronic operator V. In terms 
of the vibrational density of states 

D^(e)=Y^S{e-e^) (7) 

ft 

and thermodynamic weighting factor 

P W (e)= p|)| e=4i)! (8) 

the dephasing between two chromophores given in Eq. [5] 
becomes 

ai a*(t) = a^(Q) JdeJ de'P^{e)D^(e)P^\e')D^(e')x 

e -i(Ei+e-Ej-e')t 

(9) 

while for i = j, unitarity of U^(t,t ) ensures that 
dia*(t) = a;a*(0). 
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If the vibronic term V is assumed to be purely elec- 
tronic, so that V = Y^ij VijfJtf fij, then 

Vifij-m = Vij (0| mvntufv+t |0) (0| v n /'.,/'„, /' ,' |0) JJ 5 nkl7> 

(10) 

where the matrix elements are Franck- Condon overlaps 
between vibrational states of the excited and unexcited 
chromophores. 

Further simplification can be achieved by replacing the 
oscillatory integrals in Eq. [6] with their average over some 
period At. As At grows, this average goes to 1 if the 
energies in the exponential sum to 0; otherwise. In- 
cluding the effects of Lorentzian dephasing between chro- 
mophores in this average yields an exponential linewidth, 
so that 



(11) 



At J ~~" /3 2 + t 2 

Eq. [6] then becomes 



dt 

where 



d 

— (did*) = -iV ik T k ^i{a k a*) + i{a i a* k )V k jT k ^ j (12) 



Ti 



t->3 



de I de'D^(e)D^(e') 

P (l) (e)F(e, e')h{E l + e - Ej - e'), 



(13) 



F(e, e') = (0| iw^ti |0) (0| v^u+^j |0> | e=e(i) e , =eW 

(14) 

is the Franck Condon overlap and h(AE) is the 
linewidth. 



D. The thermalized density matrix 

In Eq. [121 T^j plays the role of the overlap be- 
tween donor emission and acceptor absorption spectra in 
Forster theory: it reflects the number of acceptor states 
sufficiently close in energy to an occupied donor state to 
receive population. The asymmetry of T^j reflects irre- 
versible flow of population from high- to low energy chro- 
mophores: for Ei > Ej, every potential donor state on 
chromophore i has a degenerate acceptor state on chro- 
mophore j, while the reverse is true only if > Ei — Ej. 

The asymmetric form of Eq. [12] which results from 
this irreversible flow can be rectified by defining a sym- 



metrized density matrix pij = e ai+aj aid* and effective 
coupling Vij = VijTi^je" 1 ^"*! . Requiring that the effec- 
tive coupling be Hermitian yields ^r^- — e 2li ~ 2lj , while 
requiring that Eq. [T2] be unaffected by the change of 
variables requires that a, = —74. Performing this sym- 
metrization and approximating the dephasing between 



chromophores as exponential yields a master equation of 
the form 

^p = -i[V ) p}-(l-8 ij )Tp. (15) 

For a constant density of states and Boltzmann proba- 
bility distribution, the Lorentzian dephasing is well ap- 
proximated by setting T = ksT. 

A byproduct of this symmetrization procedure is that 
Pij now incorporates information about the equilibrium 
distribution. In Eq. Q2] with all coherence terms set to 0, 
equilibrium is given by the condition of detailed balance, 
when 



T, 



djdj T^j 



(16) 



for all i,j pairs. Due to the symmetrization procedure, 
this condition is satisfied when p is proportional to the 
identity matrix. In this way, the current approach dif- 
fers significantly from that of Haken, Reineker and Strobl 
0, llOj . which applies a quantum master equation of 
the form of Eq. [15] to an unsymmetrized density ma- 
trix, yielding equilibrium distributions in which all chro- 
mophores have equal occupations, regardless of excita- 
tion energy or ambient temperature. If the density of 
states is again assumed to be constant and (e) to be 
Boltzmann, 014 = —"fi = fiEi/2, so that at equilibrium, 

aid; 



In the high temperature limit, where the weighting fac- 
tors approach unity, the thermalized and unthermalized 
density matrices approach one another. For this reason, 
Eq. [15] may be thought of as a generalization of the 
HRS model which describes dynamics at low tempera- 
tures. Due to the identical form, many results of the 
HRS model will carry through unchanged to describe the 
dynamics of the thermalized density matrix. 

A second consequence of the symmetrization procedure 
is that thermodynamic weighting factors have been incor- 
porated into the effective coupling. Assuming a constant 
density of vibrational states with Boltzmann weighting 
factors and approximating the linewidth and Franck Con- 
don factors as ¥(e,e%(Ei+e-Ej-e') = 5(Ei+e-Ej-e') 
yields Vij = Vije~ l3 ^ Ei ~ Ej ^ 2 , so that the effective cou- 
pling strength decays exponentially with the difference 
in excitation energies at a rate dependent on the temper- 
ature. Because ksT is only 200 wavenumbers at 300K, 
the effective interaction between two chromophores will 
decrease rapidly even for small separations in excitation 
energy. Because of this, the effective coupling matrix V 
will be both weaker and effectively more sparse than the 
matrix of bare electronic couplings, V. 



II. THE TWO CHROMOPHORE SYSTEM 

The equations of motion derived in the previous section 
jive a simple explanation for the long coherence lifetimes 
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observed in [H-Q]- For a system of two chromophores, 

—AE ) ' appl y in S Ec l- El twice yields 
differential equations for density matrix components 



with V 
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Here the total population remains constant, while the 
sum of off diagonal elements decays as e~ rt , where T = 
(25fs)~ 1 at 300K. However, both the population imbal- 
ance and the difference of off diagonal elements behave 
as damped harmonic oscillators. If If P = (poo — Pu) or 
(pio - Poi), P(t) = Ae x + l + Be x -\ with 



-r ± vr 2 - i6d 2 



(17) 



In the underdamped limit when T < 4d, P(t) oscillates 
within the exponential envelope e~ rt / 2 . When V » Ad, 
the system is overdamped and decays without oscillating 
at two rates, with A+ — > — 4d 2 /T, and A_ — > — T. 

Significantly, long coherence times do not depend on 
isolation of the system from the environment or a slow 
rate of dephasing T. On the contrary, long coherence 
times are obtained in both the underdamped limit, when 
r — > and in the overdamped limit when T » 4d. 
Because of this, long coherence times will be observed in 
both the high and low temperature limits. 

The preceding analysis has been previously derived in 
the context of the HSR model d,[loj]. However, the ther- 
modynamic weighting of the coupling matrix elements 
is crucial to accurately describing the timescale for exci- 
ton transport or the decay of coherence. The exponen- 
tial weighting term Ti^je li ~' lj which multiplies the bare 
matrix element Vij, can easily change a particular chro- 
mophore pair from the under- to the overdamped limit. 
For a constant density of states and Boltzmann weight- 
ing factors, this term evaluates to e~^ Ei ~ Ej ^ 2 , so that 
chromophores which are distant in energy relative to the 
ambient temperature will have very weak effective cou- 
plings. 

A final observation about the two chromophore system 
is that extremely long coherence times may correspond 
to slow exciton transport through the antenna complex. 
Because the population imbalance (p"oo — Pu) and the 
antisymmetric off diagonal component (pio — poi) obey 
the same differential equation, long coherence times cor- 
respond to slow decay of population imbalances. As the 
differential equations decouple from one another, this is 
not an example of coherence affecting exciton transport; 
rather, coherence and population dynamics are paramet- 
rically varying with respect to the same parameters. 



III. COMPARISON WITH EXPERIMENT - 
PE545 

The principles which govern the evolution of the two 
chromophore system apply as well to the case of a pho- 
tosynthetic complex, with the exception that, with more 
than two chromophores, there are now multiple pathways 
for the exciton to follow through the complex. The effect 
of these pathways is limited, however, by the thermo- 
dynamic weighting of the effective coupling, which may 
make the effective coupling matrix both weaker and ef- 
fectively more sparse than the bare electronic coupling 
matrix. 

In order to test the theory derived in Section HI Eq. 
US] was used to describe dynamics in the PE545 antenna 
complex of the cryptophyte algae Rhodomonas CS24, 
which was found to display long coherence lifetimes in 
1- 

Because the effective coupling strength is exponen- 
tially dependent on the energy spacing between two chro- 
mophores, it is important to use accurate excitation en- 
ergies in order to yield the correct dynamics. Here the 
error in Hartree Fock or ci singles calculations, which 
may be sizeable fractions of an electron volt, may prove 
unacceptably large. For this reason, the effective cou- 
plings and site energies for PE545 were taken from [2lj j . 
where the couplings were calculated ab initio, but the 
site energies were found by matching to several static 
and dynamic spectra, as calculated using Rcdficld theory. 
The theory was then tested by comparing calculated line- 
shapes and times for excitons to transfer between chro- 
mophores to those measured in experiment. Here it must 
be noted that the procedure of finding excitation ener- 
gies by matching to excitation spectra may yield an ar- 
tificially good agreement for lineshapes calculated using 
a different theory but the same excitation energies. 

The validity of the assumptions made in deriving Eq. 
1151 - namely, the assumption of a continuous vibrational 
density of states and a delta function linewidth, can 
be tested by constructing an effective density of states, 
corresponding to a discrete vibrational spectrum with 
a finite linewidth caused by dephasing between chro- 
mophores. Figure [5] shows the effective density of states 



D ett (co) 



-fi\u-Ui(n+l/2)\ 



n.i 



(18) 



calculated at t = 77 K and 300A" for the set of vibrational 
energies taken from [2l|. At both temperatures, the ef- 
fective density of states levels out at vibrational energies 
of a few hundred wavenumbers. At both temperatures, 
the effective density of states shows a dip at zero energy, 
due to the zero point energy of the vibrational modes. 
This drop is more accentuated at low temperatures, due 
to the narrowness of the resulting linewidth. 

Because of this, the assumption of a constant density of 
states may break down for pairs of chromophores which 
are close to degenerate, or at low temperatures. In these 
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the dipole correlation function 

d(t) =< ft{t)p(o) > . 



(19) 



400 600 800 
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FIG. 2: Effective vibrational density of states, calculated for 
PE545 assuming linewidth h(AE) = e _/3|AB| . 



limits, Ti^j should be calculated using a finite summa- 
tion in place of the integrals in Eq. Q21 

The description of dynamics in PE545 was tested by 
calculating lineshapes for absorption, circular dichroism 
and fluorescence spectra, as well as the times necessary 
to transfer excitation between pairs of chromophores. 

PE545 consists of six phycoerythrobilin (PEB) and two 
dihydrobiliverdin (DBV) chromophores, held in place by 
a dimer of two a/3 monomers. Each monomer contains 
three PEB chromophores on the j3 subunit and one DBV 
chromophore on the a subunit. The DBV chromophores 
are redshifted with respect to the PEB chromophores, 
with absorption maxima at 569 nm, compared to 545 for 
the PEB chromophores [13] . Fig Q] shows the spatial ar- 
rangement of the chromophores in the complex. Excitons 
escape the comples by first making their way to the DBV 
bilins, then slowly settling onto the lowest energy DBV 
bilin before escaping the complex. Fluorescence exper- 
iments [23| indicate that excitons primarily escape the 
complex via a single DBV bilin. Due to spectral over- 
lap in the DBV and PEB bands, it may be difficult for 
spectroscopic experiments to distinguish precisely which 
chromophores are excited. 



Static Lineshapes — Absorption, Fluorescence 
and Circular Dichroism 



Spectra for absorption, circular dichroism and fluores- 
cence can be calculated from time dependent expectation 
values of a density matrix which has been supplemented 
by addition of the ground state [2J] . Here the quantities 
of interest, (ajtip(i) can be found by removing the ther- 
modynamic weighting factors from the thermalized den- 
sity matrix, in which the dynamics are calculated. The 
absorption spectrum is given by the Fourier transform of 



where the dipole operator ft = X^A^'GO ( u l + |0) (*!)• in 
terms of the density matrix, this is given by 



d(t) = Trtfftt)), 



(20) 



where p(0) = V(0)®^*(0), <M°) = 1°) + 10 and 

pit) is found by propagating p according to Eq. [T5] The 
parameter e factors out of the Fourier transform and is 
discarded upon normalization of the spectrum. Dephas- 
ing between the ground and excited states is treated as 
an exponential decay term To; = (400cm _1 ), taken from 

& 

In the same way, the circular dichroism spectrum is 
found by taking the Fourier transform of 



m(t) = < m(i)/2(0) > 



(21) 



where m — fix R is the magnetic dipole operator. Figures 
[3] and 3] compare the calculated absorption and circular 
dichroism lineshapes to those measured in experiment 
in [2l[. Because the current theory does not account 
for reorganization energy of the chromophores following 
excitation, these spectra have been shifted in energy so 
that the peak of the calculated spectrum aligns with the 
peak of the experimental spectrum. 

In contrast to the absorption and circular dichroism 
spectra, the fluorescence spectrum reflects emission at 
large times, for a system which has settled into the low- 
est energy eigenstate. If this eigenstate is given by ip e 
and ^(0) = ijj e + e |0), the fluorescence spectrum is pro- 
portional to w 3 \f(id)\ 2 , where 

f{t) =< ftp(t) > (22) 

and p(0) = "0(0) ® V'(O)- Figure [5] compares the normal- 
ized fluorescence spectrum to that measured in experi- 
ment. 

At 300K, the calculated lineshapes show good agree- 
ment. The width of the absorption and the negative lobe 
of the circular dichroism spectra agree closely with the 
experimental widths, although the positive lobe of the 
circular dichroism spectrum has a broader tail and less 
accentuated peak than the calculated lobe. The calcu- 
lated fluorescence spectrum agrees with experiment in 
the peak region and on the blue side of the maximum, 
but falls off more quickly than experiment on the red side 
of the maximum. 

Agreement with experiment is somewhat worse at 77K. 
The absorption spectrum gives the correct width in the 
peak region, but falls off more quickly than experiment 
on the red side of the maximum and less quickly than 
experiment on the blue side. The fluorescence spectrum 
shows good agreement with the experimental lineshape 
on the blue side of the maximum, but falls off much more 
quickly than the experiment on the red side. Both the 
absorption and fluorescence spectra miss a broad tail on 
the red side of the fluorescence. The circular dichroism 
spectrum at 77K was not available for comparison. 
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FIG. 3: Normalized absorption spectra uu * d(ui) for PE545, 
compared to experimental values taken from [21|. a) 300K, 
b)77K. 



FIG. 5: Normalized fluorescence spectra ui s * d(uA for PE545, 
compared to experimental values taken from [21]. a) 300K, 
b)77K. 



- Experiment 
— Theory 

Theory (unshifted) 




18000 19000 20000 21000 

Frequency (cm -1 ) 



FIG. 4: Normalized circular dichroism spectra oj * m(u) cal- 
culated at 300K compared to experimental values taken from 

El. 



B. Exciton Transfer Times 

While static lineshapes incorporate some dynamical 
information in the form of a Fourier transform of a 
time dependent correlation function, this information is 
obscured somewhat by the thermodynamic factors fa- 
voring population of low energy chromophores and by 
finite linewidths due to the decay of coherence. In 
[25| ,timescales for exciton transfer were measured more 
directly by finding Evolution Associated Difference Spec- 
tra (EADS) associated with particular transfer lifetimes 
following an excitation pulse of a particular frequency. 
Here interpretation is complicated by spectral overlap be- 
tween chromophores, and by the difficulty in distinguish- 
ing between similar transfer lifetimes, so that only a few, 
picosecond scale lifetimes could be measured, and the 
chromophores being populated/depopulated could not be 
completely identified. 

Within these limitations, the calculated rates for exci- 
ton transfer show good agreement with experiment. Ta- 
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FIG. 6: Experimentally measured transfer times, taken from 
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P&2D 

Pl5SD 7.7 

FIG. 7: Transfer / coherence decay times (ps) shorter than 
100 ps for PE545 calculated at 77K. 



ble |6] shows the experimentally measured EADS transfer 
times, measured at 300K and 77K, for excitation wave- 
lengths of 485 nm and 530 nm. Tables [8] and [7] show 
experimental decay times A+ 1 , calculated using Eq. [17] 

As with static spectra, agreement between theory and 
experiment is somewhat better at 300K than at 77K. At 
300K, population of a DBV chromophore is found to oc- 
cur in 150 fs, vs an experimentally measured 250 fs, while 
transfer of population between DBV chromophores is cal- 
culated to occur in 20.4 ps, compared to an experimental 
value of 23.4 ps or 16.4 ps, depending on the excitation 
frequency. Several pairs of chromophores are calculated 
to have transfer times in the range of 1-3 ps, comparable 
to the remaining experimental rate of 1.84 ps. 

Agreement with experiment is considerably worse at 
77K. Population of a DBV chromophore is now calcu- 
lated to occur in 956 fs, compared to an experimental 
value of 960 fs, while a 3.7 ps time for transfer between 
/3soc and /?i58C matches well with a 3 ps experimental 
transfer times. However, the 8.4 ps calculated for trans- 
fer between DBV chromophores and the 7.7 ps calculated 
between j3^D and f3%2D do not correspond well with ex- 
perimental times of 30 ps and 17.7 ps. Finally, the theo- 
retical calculations include an underdamped transfer be- 
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FIG. 8: Transfer / coherence lifetimes (ps) shorter than 100 
ps for PE545 calculated at 300K. 



tween a±gB and Ps2C which is not seen in the experiment, 
perhaps because the real part of the transfer time, at 60 
fs, is shorter than the experimental resolution of 120 fs. 

C. Density Matrix Propagation vs. Pairwise Decay 
Model 

The agreement between the calculated pairwise decay 
times and transfer times measured in experiment raise 
the question of which information in p is necessary to 
describe the dynamics of excitons in the complex. In 
principle, the entire density matrix is necessary to de- 
scribe these dynamics, as excitons may follow any path- 
way through the complex, so that the amplitude to be on 
'liny particular chromophore is the sum of many interfer- 
ing pathways. However, due to thermodynamic weight- 
ing of the effective coupling matrix, many of these path- 
ways are strongly overdamped, with weak effective cou- 
pling and slow decay of any population imbalance. Be- 
cause of this, the network of efficient transfer pathways 
is relatively sparse. This can be seen in Tables [8] and [Jj 
or pictorially in Figure |9] where the opacity of the line 
connecting two chromophores reflects the rate of decay 
for population imbalances (pa — pjj) (or antisymmetric 
coherence terms (pij — hi)) between them. As a result of 
this sparsity, there are relatively few efficient pathways 
through the complex which may interfere with each other 
effectively. In Figure [9] two interfering pathways would 
form a closed curve of lines with nearly equal opacity. 

At 77K, aigA, Psnc and 0*193 form the only closed 
curve, while none are present at 300K (the 150 fs transfer 
time between a\gs and /3&2C is much faster than the 2 
ps transfer times between aigs,/3i58C and /fee, /3i5sc)- 

Due to the lack of interfering pathways, the flow of 
excitons through the PE545 antenna complex is well ap- 
proximated by a minimal model, in which population im- 
balances in the thermalized density matrix decay expo- 
nentially as 

-tl {pa - hi ) = ~ A + (hi - hi ) > ( 23 ) 

with A+ given by Eq. [17] 

The close agreement between the minimal model and 
the full density matrix propagation is shown in Figures 
ESI EE HI E2 HI EE5J HE and E3 which compare the 
two models for initial conditions with the excitons local- 
ized on each individual chromophore. The largest depar- 
ture between the two models is observed for 0:19,4, /?82C 
and 0195 initial conditions at 77K, which correspond to 
the closed curve of efficient transfer pathways previously 
identified. 

The close agreement between the two models explains 
why the experimental transfer rates show good agree- 
ment with the pairwise transfer rates, and indicates that 
the information contained in the long lived antisymmetric 
coherence terms, which are not included in the exponen- 
tial decay model, does not play a vital role in the trans- 



9 




FIG. 9: Rates of population transfer / coherence decay in 
PE545. Lines connecting chromophores have opacity propor- 
tional to AT. 1 , with full opacity corresponding to AT. 1 = 1 ps. 
a) Rates calculated at 77K b) Rates calculated at 300K. 



for of exciton population. Although the antisymmetric 
coherence terms show the same picosecond timescales 
for decay as the population imbalances, the lack of ef- 
ficient, closed cycles in the PE545 complex leave little 
opportunity for interference. In addition, all initial con- 
ditions settle rapidly to the same equilibrium distribution 
at 300K, so that any "quantum speedup" would offer 
only transient advantages. For these reasons, it appears 
that the importance of coherence terms in PE545 is dy- 
namic rather than informational - coherence terms are 
important because they are required to give the correct 
transfer rates between chromophores, but the informa- 
tion contained in these terms does not appear to play an 
essential role in the transfer of the exciton population to 
the radiating DBV chromophore. 




Time (ps) 

FIG. 10: Exciton populations vs time calculated at a) 300K 
b) 77K for initial state aigs- Solid lines give populations 
calculated by Eq. 1151 while dashed lines give populations 
calculated using the minimal exponential decay model. 



IV. CONCLUSIONS 

This paper has presented a new theory of coherent 
dynamics in photosynthetic complexes. Strong interac- 
tion with a reservoir causes rapid thermalization of vi- 
brational states, allowing equations of motion for slowly 
varying electronic coefficients to be derived. By incor- 
porating thermodynamic quantities into the definition of 
a thermalized density matrix and an effective coupling, 
these equations of motion are reduced to the Hakcn, 
Reineker, Strobl form. Unlike the original HRS model, 
the incorporation of thermodynamic information causes 
the system to settle to the correct thermodynamic equi- 
librium rather than to equal population of all states. The 
resulting theory gives good theoretical lineshapes for ab- 
sorption, fluorescence and circular dichroism for PE545 
at 300K and 77K, and yields exciton transfer times which 
agree with those found by EADS experiments. For both 
kinds of experiments, agreement is better at 300K than 
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FIG. 11: Exciton populations vs time calculated at a) 300K 
b) 77K for initial state psoc- Solid lines give populations 
calculated by Eq. 1151 while dashed lines give populations 
calculated using the minimal exponential decay model. 



FIG. 12: Exciton populations vs time calculated at a) 300K 
b) 77K for initial state aisA- Solid lines give populations 
calculated by Eq. 1151 while dashed lines give populations 
calculated using the minimal exponential decay model. 



at 77K. 

Long lived coherences between chromophores, such as 
those observed in [l]-|ij], are simply explained in terms of 
an overdamped harmonic oscillator, in which the strength 
of the damping is proportional to the temperature and 
the effective coupling decreases exponentially with the 
excitation energy separation between two chromophores. 
Because of this, coherence terms may survive for arbi- 
trarily long times, even in the high temperature limit for 
a system which interacts strongly with its environment. 

Somewhat surprisingly, the long lifetime for survival 
of coherence terms appears to have a minimal effect on 
the transfer of exciton population in the PE545 antenna 
complex. This can be understood as resulting from the 
relatively sparse network of efficient transfer pathways 
in this complex, which has the effect of limiting inter- 
ference between competing pathways. A minimal expo- 
nential decay model describes the flow of exciton popula- 
tion through this network accurately at 300K, when this 



network has no closed cycles, but departs from the full 
density matrix calculation at 77K, when such a cycle al- 
lows for interference between multiple pathways through 
the complex. It is possible that other photosynthetic 
molecules, with functions more sophisticated than the 
simple direction of exciton population to the lowest en- 
ergy chromophore, may make more sophisticated use of 
this information. 

The method of thermalized reduction of the density 
matrix introduced in this paper allows the derivation of 
quantum coherent equations of motion in a regime of high 
temperature and strong interaction with the surrounding 
environment which is often considered inimical to quan- 
tum mechanical behavior. Rather than being destroyed 
rapidly by interaction with the reservoir, coherence terms 
may persist or even be preserved by interaction with the 
reservoir. The conditions necessary for reduction - that 
the reduced degrees of freedom thermalize rapidly with 
respect to the evolution of the unreduced degrees of free- 
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dom - are satisfied in the limit of strong interaction with 




the reservoir, while the thermodynamic weighting of the 
effective coupling matrix may yield slow evolution of the 
unreduced degrees of freedom even when the bare cou- 
pling is relatively strong. The simplicity of the conditions 
necessary for the thermalized reduction procedure show 
that quantum mechanical behavior does not require low 
temperatures or elaborate isolation from the surround- 
ings to be observed. 
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FIG. 16: Exciton populations vs time calculated at a) 300K 
b) 77K for initial state P%2d- Solid lines give populations 
calculated by Eq. 1151 while dashed lines give populations 
calculated using the minimal exponential decay model. 
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FIG. 17: Exciton populations vs time calculated at a) 300K 
b) 77K for initial state /3i58D- Solid lines give populations 
calculated by Eq. 1151 while dashed lines give populations 
calculated using the minimal exponential decay model. 
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